PROCEEDINGS 

OF SCIENCE 



MKPH-T-07-15 



Relativistic corrections to the static potential at 

0(1 /m) and 0(l/m 2 ) 



Yoshiaki Koma 

Numazu College of Technology, Numazu 410-8501, Japan 



E-mail: koma@numazu-ct .ac.jp 



Miho Koma 1 

Institut filr Kernphysik, Johannes Gutenberg-Universitat Mainz, D-5 5099 Mainz, Germany 
E-mail: }nkoma@kph . uni-mainz . de 



Hartmut Wittig 

Institut filr Kernphysik, Johannes Gutenberg-Universitat Mainz, D-5 5099 Mainz, Germany 



E-mail: i*ittig@kph . uni-mainz . de 



We investigate the relativistic corrections to the static potential, i.e. the 0(1 /m) potential and 
the 0(l/m 2 ) velocity-dependent potentials, in SU(3) lattice gauge theory. They are important 
ingredients of potential nonrelativistic QCD for heavy quarkonium. Utilizing the multi-level al- 
gorithm, we obtain remarkably clean signals of these potentials up to r = 0.9 fm. We observe 
long range nonperturbative contributions to these corrections. 



The XXV International Symposium on Lattice Field Theory 
July 30 - August 4 2007 
Regensburg, Germany 



*Speaker of "Determination of the velocity-dependent potentials at 0(l/m 2 )" 
' Speaker of "Relativistic correction to the static potential at 0( 1 //n)" 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http: //pos . sissa . it/ 



Relativistic corrections to the static potential at 0(1 jrri) and 0(1 /m 2 ) 



1. Introduction 

A possible strategy of studying heavy quarkonium in QCD is to employ potential nonrelativis- 
tic QCD (pNRQCD) |jl|, g[ ||, Q], which is obtained by integrating out the scale above the heavy 
quark mass m 3> Aqcd and the scale mv, where v is quark velocity. The effective hamiltonian of 
pNRQCD up to 0(1 /m 2 ) [g] is then given by 

H = lL + lL +V (0)( r) + _L v (i,o) (r) + J_ v (o,i) (r) 
2m i 2rri2 m\ ni2 

+ _Ly(2,0) (r) + J_y(0,2) (r) + _J_y(l,l) (r) + 0(1/m 3 ) ; (U) 
m\ m-2 m\m2 

where m\ and ni2 denote the masses of quark and antiquark, placed at f\ and F2, respectively. 
The static inter-quark potential V^(r = \r\ — ~?2 | ) emerges, accompanied by relativistic corrections 
classified in powers of l/m. The potentials V^ L0 \r) = V^ l) (r) (= V^\r)) are the corrections at 
0(1 /m). The potentials V^ 2 '°\r), V {0 - 2 \r), and V^(r) are the corrections at 0(1 /m 2 ), which 
contain the leading order spin-dependent potentials [§, 0] and the velocity-dependent poten- 
tials [|| |9p. Spin-dependent potentials are relevant to describing the fine and hyper-fine splitting 
of heavy quarkonium spectra. Once these potentials are determined, various properties of heavy 
quarkonium, not only the ground state but also excited states, e.g. full spectrum and wave functions, 
can be investigated systematically by solving the Schrodinger equation. 

One may rely on perturbation theory to determine these potentials in the short-distance region. 
For instance, perturbative studies of the 0(1 /m) potential yield V^ l '(r) = —CFCAOC 2 /(4r 2 ) [JTO , U 



§], where C F =4/3 and C A = 3 are the Casimir charges of the fundamental and adjoint representa- 
tions, respectively, and a s = g 2 / (4n) the strong coupling (for the expression beyond leading-order 
perturbation theory, see [|2|]). However, since the binding energy of a quark and an antiquark is 
typically of the scale mv 2 , which can be of the same order as Aqcd due to the nonrelativistic na- 
ture of the system, v <C 1, as well as the fact that perturbation theory cannot incorporate quark 
confinement, it is essential to determine the potentials nonperturbatively. 

Monte Carlo simulations of lattice QCD offer a powerful tool for the nonperturbative determi- 
nation of the potentials. The static potential can easily be determined from the expectation value 
of the r x t w rectangular Wilson loop by V^°\r) = — lim tw ^oo(l / t w )ln(W (r,t w )) . The result is well 
parametrized by the Coulomb plus linear confining term, 

V^(r) = -- + ar + n, (1.2) 
r 

where c denotes the Coulombic coefficient, a the string tension and /I a constant term. Recently, 
we investigated the 0(1 /m) potential [13] and the 0(1 /m 2 ) spin-dependent potentials [|l4|, 15] 



on the lattice utilizing a new method, and obtained remarkably clean signals up to distances of 
around 0.6 fm. We then observed a certain deviation from the perturbative potentials at intermediate 
distances. 

In this presentation we further investigate the relativistic corrections to the static potential. In 
particular, we aim to clarify the long-distance behavior of the 0(1 /m) potential, and to determine 
the 0(1 1 m 2 ) velocity-dependent potentials. One attempt to determine the velocity dependent po- 
tentials on the lattice was published a decade ago [|16|], but the signal was lost due to large statistical 
errors. 
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2. Spectral representation of the 0(1 /m) and 0(1 /m 2 ) potentials 

According to pNRQCD, the 0(1 /m) and 0(1 /m 2 ) potentials can generally be expressed by the 
matrix elements and the energy gaps of the spectral representation of the field strength correlators 
on the quark-antiquark source [§, [3j. 

Writing the eigenstate of the pNRQCD hamiltonian at 0(m Q ) in the 3 <g>3* representation of 
SU(3) color, which corresponds to the static quark-antiquark state, as \n) = \n;r\,ri), the correlator 
of two color-electric field strength operators E' = (i = 1,2,3), put on f a and ft, (a,b = 1,2) in 
space and separated t = t\ — ?2 in time, takes the form 

CO 

C(r,t) = £(0|^ ! '(r fl )|«}(n|^(r,)|0)^( A£ "»W , (2.1) 

n=\ 

where AE nQ (r) = E n (r) —Eo(r) denotes the energy gap with E (r) = V^(r). 
Then, the 0(1 jm) potential is given by 

W- 2 VL {AEno{r)) 2 . M 

where two color-electric field strengths are attached to one of the quarks. 
The spin-independent part of the 0(1 /m 2 ) potential is written as 

m\\L p r z ' 



+ ~2 «> + ^ (r)h 2 + (r) 

m\\l p r L 1 

1 f 1 ri.i>,„ 1 



m\m2 



^{/^•K,v; ( ^ 1) (r)}-^ M) (r)(r 1 .r 2 +r 2 .r 1 )+v r (M) (r)) , (2.3) 



where l a = r x p a . The functions specified by the subscripts p 2 and I 2 are related to the velocity- 
dependent potentials, Vb, V c , Vd and V e defined in Refs. by 



V£°\r) = Vp { r) = V d (r) - \v e (r) , V™(r) = V™(r) = V e (r) , 

Vy\r) = -V b (r) + 2 -V c (r), V^\r) = -V c (r) , (2.4) 



and then we have 



w-|«,£ (°I»*WMW J MP) , (2 . 5) 

„.(,) _ 3 (^4) £ W'WyWW , (2.6, 



Ur) = 3*' 5 (A£^oP ' (2 - 7) 



3/r,r ; - $A f, (0\gE\n)\n)(n\gE](n)\0) 
Ve(r) = ~2\r r ~^) h • (2 - 8) 
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For Vd and V e two field strengths are attached to one of the quarks as for the V^, while for Vb and 
V c two field strengths are attached to the quark and the antiquark, respectively. 

Thus, once the matrix elements and the energy gaps are known from the behavior of the field 
strength correlators on the quark-antiquark source, one can compute the potentials. 



3. Numerical procedures 

We work in Euclidean space in four dimensions on a hypercubic lattice with lattice volume 
V = L?T and lattice spacing a, where we impose periodic boundary conditions in all directions. 
We use the Polyakov loop correlation function (PLCF, a pair of Polyakov loops P separated by a 
distance r) as the quark-antiquark source and evaluate the color-electric field strength correlators 
on the PLCF, 

C(r,t) = ((gE^ra^gE^h^c = ((gE^r^gE 1 (7 b ,h)))- ((gE 1 (r a ))) ((gE^(r b ))) , (3.1) 



using the multi-level algorithm [13, 15], where the double bracket represents the ratio of expec- 



tation values ((•••)) = (•■•)«»/ (PP*(r)}, while (■■■} PP , means that the color-electric field is con- 



nected to the Polyakov loop in a gauge invariant way. The subtracted term on the r.h.s. of Eq. ( j3.1| ) 
can be nonzero as the color-electric field is even under CP transformations. The spectral represen- 
tation of Eq. ( |3~l| ) derived with transfer matrix theory reads [15] 



OO rji 

C(r,t)=2£(0\gE\r a )\n)(n\gE\r h )\0)e-^ T / 2 ^ , (3.2) 

n=\ 1 

where the last term represents terms of exponential factors equal to or smaller than exp(— (AE\q)T), 
which can be neglected for a reasonably large T. Thus, once Eq. ([OJ) is evaluated via Monte Carlo 
simulations, we can fit the matrix element (0\gE l (r a )\n) (n\gE-i (r b )\0) and the energy gap AE„o in 



Eq. (3.2), both of which are finally inserted into Eq. (fL2j), etc. It is clear that Eq. (3.2) is reduced 



to the form like Eq. (2. 1) in the infinite volume limit T — > °°. 

We define the lattice color-electric field operator, ga 2 E'(s), from the traceless part of [U4i(s] — 
Ul t (s)]/(2i) with two-leaf modification (an average of /^(s) and F^j — i)), where U^ v (s) is a 
plaquette variable defined on the site s. We multiply the Huntley-Michael factor $17\\ on the PLCF, 
Ze [JTJ], to the lattice color-electric field to cancel the self energies at least at 0(g 2 ). 

We point out several advantages of our procedure which enable us to reduce numerical errors. 
In earlier studies of the relativistic corrections on the lattice, the Wilson loop has been used as 
the quark-antiquark source, since the corrections have been expressed as the integral of the field 
strength correlators on the Wilson loop with respect to the relative temporal distance between two 
field strength operators, t (see, e.g. [Q|). They are, in principle, measurable on the lattice, and the 



result is reduced to Eq. ( |2.2| ) once the spectral decomposition is applied by using transfer matrix 
theory, and the temporal size of the Wilson loop is taken to infinity. In practice, however, the inte- 
gration of the field strength correlator and the extrapolation of the temporal size of the Wilson loop 
to infinity cause systematic errors. In contrast, we can avoid these systematic errors once the field 
strength correlator is determined accurately, since we directly evaluate the matrix elements and the 
energy gaps. The use of the PLCF as the quark-antiquark source and its spectral representation 
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allows us to take into account the finite-J effect automatically in the fit. We should note however 
that if one uses the commonly employed simulation algorithms, it is almost impossible to evaluate 
the field strength correlators on the PLCF, or the PLCF itself, at intermediate distances with reason- 
able computational effort, since the expectation value of the PLCF at zero temperature is smaller 
by several orders of magnitude than that of the Wilson loop. However this problem is solved by 
employing the multi-level algorithm [|^]. For details of our implementation, see [|l5|]. 



4. Numerical results 



We carry out simulations using the standard Wilson gauge action in SU(3) lattice gauge theory. 
We summarize our simulation parameters in Table [T[ The lattice spacing a is set from the Sommer 
scale ro = 0.5 fm. For a reference we compute the static potential and the force from the PLCF, 

vW(r) = -Ib.<IV W > + 0(.-<«^, j_^W = V W W-V W (r- a ) | (41) 

which are shown in Fig. 0. The fitting of the potential data at j3 = 6.00 to the functional form in 
Eq. (O) yields c = 0.2808(5), oa 2 = 0.0468(1) and \La = 0.7301(4) with X 2 /N Ai = 3.5. We note 



that the large value of X 2 /Ndi just reflects the fact that the Coulombic coefficient is not strictly 
constant as a function of r, which will be clear once the second derivative of the potential is inves- 



tigated [18, 15 1, but we ignore this effect for the moment. 



Table 1: Simulation parameters used in this study. N t $\ is the number of time slices in a sublattice and Mupd 
the number of internal update within a sublattice, both are parameters for the multi-level algorithm. 
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Figure 1: Static potential V^°\r) and the force dV^°\r)/dr as a function of r, which is improved with 
treelevel perturbation theory. The potentials of different /3 values are normalized at r = 0.5 fm. 

In Fig. ^, we show the typical behavior of the longitudinal and the transverse components of 
the color-electric field strength correlators, where the quark-antiquark axis has been taken along the 
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Figure 2: Color-electric field strength correlators at /3 = 5.85 on the 18 3 24 lattice for r/a = 5. 
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Figure 3: The 0(1 /m) potential and the force. 
x axis. We select the data at r/a = 5, j8 = 5.85 as an example. Statistical errors of the correlators are 



small enough to determine the matrix elements and the energy gaps with the fit based on Eq. (3.2). 



It is evident that Eq. (3.2) describes the data very well. Here, since it is impossible to determine 
the matrix elements and the energy gaps for all n > 1 with the limited data points, we truncated 



expansion in Eq. (3.2) at a certain n = « max - The validity of the truncation was monitored by looking 
at the reduced % 2 defined with the full covariance matrix. In all cases we found that « max = 3 gave 
the best result. 

We then computed the 0(1 /m) potential and the 0(1 /m 2 ) velocity-dependent potentials, which 
are plotted in Figs. [| and |], respectively. Note that each potential contains a constant contribution, 
which depends on p. Here we normalized the potentials at r = 0.5 fm by assuming perfect scal- 
ing behavior. This assumption is justified at intermediate distance, where the data at different 
lattice spacings fall into a smooth curve. Small discrepancies at short distances, however, can- 
not be avoided in the present simulation. The data at short distances are sensitive to the way of 
discretization and the definition of the renormalized color-electric field operator. 
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Figure 4: The 0(l/m 2 ) velocity-dependent potentials 
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The functional form of the 0(1 /m) potential is not yet established nonperturbatively. In our 



previous study Q131, we found that a 1/r function describes the data well, but this result was valid 
only up to r = 0.6 fm. Now we have further long distance data up to distances of r = 0.9 fm. We 
empirically examined various functional forms and found that, if we include the data r > 0.6 fm, 
the 1/r function is not supported by the fit, while the perturbative 1/r 2 function with the linear 
term can fit the data well. Using the functional form 



tit 



+ o'r + ii' 



(4.2) 



we obtain c' = 0.090(5), a'a 3 = 0.0024(1), and \i'a 2 = 0.389(1) with £ 2 /7Vdf = 0.68 for the data 
at j3 = 6.00, where the distance r/a = 2 is omitted. In Fig. |3[ we also plot the force of the 0(1 /m) 
potential defined by 



rfy("(r) 



dr 



tion — - 



V) 



dr 



(y(i) ( r ) _y(i) ( r _ a )) / a . We find that the fit of the force to the func- 

^ + a' gives c' = 0.095(5), a' = 0.0024(1) with x 1 /N Af = 0.62, which is consistent 
with the fit result of the potential. 

By taking into account the masses (l/m\ + 1 /m-i) in Eq. (|1.1[), we may estimate the correction 
to the string tension in the static potential. For charmonium, m c = 1.3 GeV, we find (2/m c )o' = 
0.179 GeV/fm, which is compared to a = 1.07(1) GeV/fm. The correction is about 17 %. For bot- 
tomonium, nib = 4.7 GeV, we find (2/nib)o' = 0.049 GeV/fm, so that the correction is about 5 %. 

Next, we may characterize the functional form of the velocity-dependent potentials. Here, 
motivated by the minimal area law (MAL) model [M S], we fit the potentials to the same functional 
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Table 2: Fit result of the velocity dependent potential to the function V& t (r) = —c/r+ar + ji. 
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Figure 5: Test of the BBMP relations in Eq. (4.3). 



form as the static potential. Again, we omit the point at rja = 2 from the analysis of the data taken 
at p = 6.00. The fitting results are summarized in Table |2[ The global structure of the data seems 
to be well described by the fitting function as seen in Fig. fij though % 2 /N& may be relatively large. 

Finally we examine some nonperturbative relations which connect the velocity-dependent po- 
tentials to the static potential [ISL 0], which are often called the BBMP relations, 



, ftW+ 2V jW = -Iv<») (r) + ^, W+a , W = -£^. (43, 

These relations are derived by exploiting the exact Poincare invariance of the field strength correla- 
tor, and hence one expects corrections due to lattice artifacts, which would vanish in the continuum 



limit. These relations can be regarded as an extension of the Gromes relation [19] for the 0(1 /m 2 ) 
spin-dependent potentials. In Fig. ||, we show the result, where the constant contribution of the 
potentials is normalized at r = 0.5 fm. The BBMP relations seem to be satisfied, though we see a 
small discrepancy especially at short distances. 



5. Summary 

We have investigated the relativistic corrections to the static potential, the 0(1 /m) potential 
and the 0(1 /m 2 ) velocity-dependent potentials, in SU(3) lattice gauge theory. They are important 
ingredients of the pNRQCD hamiltonian for heavy quarkonium. 

By evaluating the color-electric field strength correlator on the PLCF with the multi-level 
algorithm, and exploiting the spectral representation of the correlator, we have obtained a very 
clean signal for these potentials up to r = 0.9 fm. The 0(1 /m) potential contains a linearly rising 
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nonperturbative contribution. The 0(1 /m 2 ) velocity-dependent potentials are non-vanishing at 
long distances. 

All potentials at different /3 values, normalized at r = 0.5 fm, show a reasonable scaling be- 
havior. The BBMP relations are apparently satisfied around r ~ 0.5 fm. Although we have applied 
the Huntley-Michael prescription in the present study to remove the self-energy contributions of 
the field strength operator, a more systematic, non-perturbative renormalization procedure of field 
strength operators is highly desirable. Since the statistical errors of the potentials are reduced sig- 
nificantly owing to the multi-level algorithm, we now face such a delicate problem. 

The comparison with various models [(9, 20] and phenomenology [^1], ^2, 23] are of course to 
be done. In particular, it is quite interesting to examine the effect of the 0(1 /m) potential on the 
spectrum as this is the leading-order relativistic correction in the Ijm expansion. 
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